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Abstract. A collection of spherical obstacles in the ball in Eu- 
clidean space is said to be avoidable for Brownian motion if there 
is a positive probability that Brownian motion diffusing from some 
point in the ball will avoid all the obstacles and reach the bound- 
ary of the ball. The centres of the spherical obstacles are generated 
according to a Poisson point process while the radius of an obsta- 
cle is a deterministic function depending only on the distance from 
the obstacle's centre to the centre of the ball. Lundh has given the 
name percolation diffusion to this process if avoidable configura- 
tions are generated with positive probability. An integral condition 
for percolation diffusion is derived in terms of the intensity of the 
Poisson point process and the function that determines the radii 
p/j • of the obstacles. 

Oh 

•s 

1. Introduction 

Lundh proposed in [10] a percolation model in the ball B = {x G 
M. d : | a; | < 1}, d > 3, involving diffusion through a random collection 
of spherical obstacles. In Lundh's formulation, the radius of an ob- 
O ■ stacle is proportional to the distance from its centre to the boundary 

§ = {x € M. d : \x\ = 1} of the ball. The centres of the obstacles 
are generated at random by a Poisson point process with a spherically 
symmetric intensity \x. Lundh called a random collection of obstacles is 
avoidable if Brownian motion diffusing from a point in the ball B has a 
positive probability of reaching the outer boundary § without first hit- 
ting any of the obstacles. Lundh set himself the task of characterising 
those Poisson intensities /i which would generate an avoidable collec- 
tion of obstacles with positive probability, and named this phenomenon 
percolation diffusion. Our main objective herein is to extend Lundh's 
work by removing some of his assumptions on the Poisson intensities 
and on the radii of the obstacles. 

Deterministic configurations of obstacles in two dimensions are con- 
sidered in detail by Akeroyd [3] and by Ortega-Cerda and Seip [32], 
while O' Donovan [11] and Gardiner and Ghergu [8] consider config- 
urations in higher dimensions. The result below is taken from these 
articles. First some notation is needed. Let B(x,r) and S(x,r) stand 
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for the Euclidean ball and sphere, respectively, with centre x and radius 
r and let B(x,r) stand for the closed ball with this centre and radius. 
Let A be a countable set of points in the ball B which is regularly 
spaced in that it has the following properties 

(a) there is a positive e such that if A, A' G A, A 7^ A' and |A| > |A'| 
then 

|A-A'| >e(l-|A|). (1.1) 

(b) there is an r < 1 such that 

B= (jB(A,r(l-|A|)). (1-2) 

AeA 

Let <p : [0, 1) — > [0, 1) be a decreasing function such that the closed 
balls {-B(A, 0(| A|)}, A G A, are disjoint, and set 

0= \Jb(\,</>(\\\). 

AeA 

Avoidability of the collection of spherical obstacles O is equivalent to 
the harmonic measure condition u(x, S, fl) > 0, where Q — B \ O and 
x is some (any) point in the domain Q. 

Theorem A. The collection of spherical obstacles O is avoidable if 
and only if 

(H < 00 if d = 2, (1.3) 



*)log((l-t)/0(t)) 
1 <p(t) d - 2 



;i - 



< 00 if d>3. (1.4) 



Our goal is is to obtain a counterpart of this result for a random con- 
figuration of obstacles. We work with a Poisson random point process 
on the Borel subsets of the ball B with mean measure dfi(x) = u(x) dx 
which is absolutely continuous relative to Lebesgue measure. (Ito 
presents a complete, concise treatment of this topic in Section 1.9 of 
his book [9]). The radius function and the intensity function v are 
assumed to satisfy, for some C > 1, 

i, V {x) < v(y) < Cu(x) 9 V 2 / 



It is also assumed that 

0(1*1) 



< c < 1 for s G B. (1.6) 



1 — |a;| 

and that 

1 - \x\)(j){x) d - 2 v{x) = () I ~ ) as I ./- 1 -> r. (1.7) 
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Let V be a realisation of points from this Poisson random point process 
and let 

A P ={jB(p,<p(p)), Q V = M\A V , (1.8) 
pev 

so that Q-p is an open, though not necessarily connected, subset of 
B. The archipelago of spherical obstacles Ap is said to be avoidable 
if there is a positive probability that Brownian motion diffusing from 
some point in Qp reaches the unit sphere § before hitting the obstacles 
Ap, that is if the harmonic measure of Ap relative to flp satisfies 
u(x, Ap, Qp) < 1 for some x in Qp. If Qp is connected then, by the 
maximum principle, this condition does not depend on x E flp. We do 
not insist, however, on the configuration being avoidable for Brownian 
motion diffusing from the origin. 

We have percolation diffusion if there is a positive probability that 
the realisation of points from the Poisson random point process results 
in an avoidable configuration. Our main result is 

Theorem 1. Suppose that (II. 5ft . (11.61) and (II .7ft hold. Percolation 
diffusion occurs if and only if there is a set of points r of positive 
measure on the sphere such that 

f- <f){x) d ~ 2 v{x) dx < oo. (1.9) 

Thus the random archipelago Ap is avoidable with positive probability if 
and only if the Poisson balayage of the measure (1 — \x\ 2 )(j)(x) d ~ 2 v(x) dx 
is bounded on a set of positive measure on the boundary of the unit ball. 

Furthermore, in the case of percolation diffusion the random archi- 
pelago Ap is avoidable with probability one. 

In the radial case the following corollary follows directly from The- 
orem m 

Corollary 1. Suppose that, in addition to (11.51) . (11.61) and (11.71) . the 

intensity v and the radius function <fi are radial in that they depend 
only on \x\. Then percolation diffusion occurs if and only if 

[ (1 - t) <j){t) d - 2 v{t) dt < oo. (1.10) 
Jo 

Lundh's result [TQl Theorem 3.1] is the case 4>{t) = c(l — t) of this 
corollary, in which case (ll.lOp becomes 

[ {l-t) d - l v{t)dt < oo. (1.11) 
Jo 

This corresponds to the condition stated by Lundh that the radial in- 
tensity function should be integrable on (0, oo) when allowance is made 
for the fact that he works in the hyperbolic unit ball. As pointed out in 
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Lundh's deduction from ( II. lip (see [TUJ Remark 3.2]) that perco- 
lation diffusion can only occur when the expected number of obstacles 
in a configuration is finite isn't correct. In fact, (11. lip holds in the case 
v(t) = (1 — t) 1_d and we have percolation diffusion. At the same time, 
the expected number of obstacles N(M) in the ball is 

E[N(M)} = [ dfi{x) = [ n ^ = oo. (1.12) 
Jb Jb I 1 — \ x \) 

Lundh's remark erroneously undervalues his work since it gives the 
impression that, in his original setting, percolation diffusion can only 
occur if the number of obstacles in a configuration is finite almost 
surely. 

The intensity u(t) = 1/(1 — t) d corresponds, in principle, to a reg- 
ularly spaced collection of points since the expected number of points 
in a Whitney cube Q of sidelength £(Q) and centre c(Q) is, in the case 
of this intensity, 

E[N(Q)\ = J dfi(x) ~ v(c(Q))Vol(Q) = ^Z^Q^d ~ constant - 

We note that there is agreement in principle between the integral condi- 
tion fll.4p for the deterministic setting and the integral condition fll.lOp 
with v{t) = 1/(1 — t) d for the random setting. 



2. AVOIDABILITY, MINIMAL THINNESS AND A WlENER-TYPE 

CRITERION 

Avoidability of a realised configuration of obstacles A-p may be rein- 
terpreted in terms of minimal thinness of Ap at points on the boundary 
of the unit ball (see [I] for a thorough account of minimal thinness). 
This is Lundh's original approach, and is also the approach adopted by 
the authors of HHiniE]. 

For a positive superharmonic function u on B and a closed subset A 
of B, the reduced function R™ is defined by 

= inf {v. v is positive and superharmonic on B and v > u on A}. 

The set A is minimally thin at r € S if there is an x in B at which 
the reduced function of the Poisson kernel P(-,r) for B with pole at 
r satisfies Rp, T Jx) < P(x,t). Minimal thinness in this context has 
been characterised in terms of capacity by Essen [7] in dimension 2 
and by Aikawa [1] in higher dimensions. Let {Qfc}^ =1 be a Whitney 
decomposition of the ball B into cubes so that, in particular, 

diam(Q fc ) < dist(Qj;,§) < 4diam(Q fc ). 

Let £(Qk) be the sidelength of Qk- Let cap (E) denote the Newtonian 
capacity of a Borel set E. Aikawa' s criterion for minimal thinness of 
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A at a boundary point r of B is that the series VF(A, r) is convergent, 
where 

Pfc being the distance from Q k to the boundary point r. A proof of the 
following proposition can be found in [SJ Page 323]. The proof goes 
through with only very minor modifications even though we do not 
insist on evaluating harmonic measure at the origin and the open set 
M\A may not be connected. 

Lemma 1. Let A be a closed subset o/B. Let 

Ai = {t G § : A is minimally thin at r}. (2-2) 

Then A is avoidable if and only if Ai has positive measure on §, that 
is if and only ifW(A,r) < oo for a set of t of positive measure on §. 

The question of whether a given set A is avoidable for Brownian 
motion is thereby reduced to an estimation of capacity. 

The following zero-one law simplifies the subsequent analysis, and 
will imply that the random archipelago is avoidable with probability 
zero or probability one, as stated in Theorem 1. Again, r is used to 
denote points on the sphere S and Ap denotes an archipelago con- 
structed as in (11.81) from a random realisation V of points taken from 
the Poisson point process. 

Lemma 2. The event that A-p is minimally thin at r has probability 
or 1. 

Proof. Whether or not the set A-p is minimally thin at r depends on 
the convergence of the series W(Ap,r). Partition the cubes {<3fc}i° 
into finitely many disjoint groups {Q k }T=i> « = 1, 2, . . ., n, so that any 
ball in Ap can meet at most one cube in each group. Then break the 
summation W(Ap, r) into corresponding summations 

\2 



W\A V , r) = Y,X l k where X k = cap (A n Q k ) (2.3) 

fc=l Pk\ T ) 

The random variables X l k in each resulting summation are independent. 
The event W*(Ap, r) < oo belongs to the tail field of the corresponding 
A£'s, hence this event has probability or 1. It follows that the event 
W(Ap, r) < oo has probability or 1. □ 

3. The expected value of the Wiener- type criterion and 

the Poisson balayage 



The proof of Theorem [T] follows the outline of Lundh's argument [10J 
and the second author's thesis [T2] . 

We work with a Poisson point process in the ball. Each realisation 
V of this process gives rise to an archipelago Ap via (11.81) . which is 
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avoidable for Brownian motion if and only if the associated Wiener- 
type series W(Ap, r) is finite for a set of r of positive measure on the 
sphere S. For a fixed r on the sphere S, the series W(Ap>, r) is a random 
variable. Proposition [T] states that its expected value is comparable to 
the Poisson balayage (jl.9|) . We denote by c and C any positive finite 
numbers whose values depend only on dimension and are immaterial 
to the main argument. 

Proposition 1. Fix a point t on the sphere §. Then 

E[W(A v ,t)]~ [ ( , ~ |X| ,] 2 0(x) d ~ 2 v(x) dx. (3.1) 



The proof of Proposition [T] depends on a two-sided estimate for the 
expected value of the capacity of the intersection of a Whitney cube 
Qk with the set of obstacles Ap in terms of the mean measure fi(Qk) 
of the cube and a typical value of the radius function <p on the cube. 

Lemma 3. For a Whitney cube Q and any point x e Q, 

E[cap {A P n Q)] ~ <P(x) d ~ 2 n(Q). (3.2) 

Lundh did not require an estimate of this type as the size of one of his 
obstacles was comparable to the size of the Whitney cube containing 
its centre. The capacity of Ap PI Q therefore depended only on the 
probability of whether of not the cube Q contained a point from the 
Poisson point process. We first deduce Proposition [1] from Lemma [3] 
and then prove Lemma [31 

Proof of Proposition^ The upper bound for E[cap (ApCiQ)} in Lemma[3] 
leads to an upper bound for the expected value of Aikawa's series ( 12. ip 
with A = Ap as follows: 



E[W(A v ,t)} = E 



\2 



, Pk(rY 

< g e^t€^" 2 mq,) 

where Xk is any point in Qk- Since the radius function is approxi- 
mately constant on each Whitney cube by (11.51) . it follows that 

E[W(A V , r)]<CJ2 [ { ] - f vt?) dx 

f (1 - \ r \ 2 ) 2 
= c / ^ l -LL- (j)(x) d ~ 2 u(x) dx. 
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In the other direction, first choose a point x k in each Whitney cube 
Q k . Then, 

/ ( 1 ~ N ,] 2 m d - 2 u(x) dx<cj2 ^ d ~ 2 ^) 

Jn \ T — x \ l Pk{T) 

<cJ2 £ -^n™v(A P nQ k )] 

k Pk ^ T > 
= CE[W(A v ,t)}, 

where the second inequality comes from the lower bound for E[cap (Apfl 
Qk)] i n Lemma [31 □ 

Proof of Lemma 0. The assumption (11.61) implies that if an obstacle 
meets a Whitney cube Q then its centre can lie in at most some fixed 
number iV of Whitney cubes neighbouring the cube Q. We label these 
cubes Q l , where the index i varies from 1 to at most N, and write Q' 
for their union. Both the distance to the boundary, and the distance to 
a specific boundary point, are comparable in Q and in Q '. Analogously, 
an obstacle with a centre in a specified cube can intersect at most some 
fixed number of neighbouring cubes. 

Consider a random realisation of points V and a Whitney cube Q k . 
By (11.51) the radius function <f> is roughly constant on the cubes Q\, say 
(f)(x) ~ 4>(xk), x G Q' k , where x k is any point chosen in Q k . Therefore, 
by the subadditivity property of capacity, 

cap^nQ,,) < C<p(x k ) d ~ 2 N(Q' k ), 

where N(Q' k ) is the number of centres from the realised point process 
V that lie in the union of cubes Q' k . Taking the expectation leads to 

E[cap (A v n Q k )] < C<j)(x k ) d - 2 E [N(Q' k )] = C <P(x k ) d - 2 fi(Q' k ). 

By (11.51) . fJ.(Q' k ) < C/j,(Qk) and the upper bound for E[cap (A-p D Qk)] 
in Lemma [3] follows. 

In the other direction we proceed, as did Gardiner and Ghergu [S], 
by employing the following super-additivity property of capacity due 
to Aikawa and Borichev [2]. Let be the volume of the unit ball. Let 
F = |J B{yk, Pk) be a union of balls which lie inside some ball of unit 
radius. Suppose also that p k < 1/ for each k and that the larger 
balls B(y k ,a~ [ 1 ^ d p k ~ 2 ^ d ) are disjoint. Then 

cap(F) >cJ2 ^P (B(y k , p k )) =c^pf- 2 . (3.3) 

k k 

Let 0o be the minimum of 4>(x) for x in Q. By ( II. 5p . 0o is comparable 
to (j){x) for any x in Q. We only consider obstacles with centres in Q 
and suppose that all such obstacles have radius <po, since in so doing 
the capacity of A-p PI Q decreases. Set 

a = mm{(£(Q)Vd)~\ {^2 d cj)^ 1 } 
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and set 

N =lA-H{Q)a l J d a 2 ' d <pT~ 1 \- 

By (11. 6p . we have a > c/£(Q). The cube Q is divided into N d smaller 
cubes each of sidelength £(Q)/N: we write Q 1 for a typical sub-cube. 
Inside each cube Q' consider a smaller concentric cube Q" of sidelength 
£{Q) I (4iV). If a cube Q" happens to contain points from the realisation 
V of the random point process, we choose one such point. This results 
in points Ai, A2, • • ., A TO , say, where m < N d . By the choice of a and N, 
each ball B(Xk,<fio) is contained within the sub-cube Q' that contains 
its centre. We set 

m 

A P , Q = \jB(X k ,(/)o). 

k=l 

Since A-p t Q C A-p H if follows from monotonicity of capacity that 
cap [A-p n Q) > cap (Ap,q) . 

To estimate the capacity of Ap q, we scale the cube Q by a. By 
the choice of a, the cube aQ lies inside a ball of unit radius and the 
radius of each scaled ball from A V q satisfies a0 o < (cd2 d )~ 1//2 . The 
only condition that remains to be checked before applying Borichev 
and Aikawa's estimate ( 13. 3 j) to the union of balls aAp^ is that the 

balls with centre a\k and radius o"J 1//d (a0o) 1-2 ^ are disjoint. They 
are if 

since the centres of the balls are at least a distance a £(Q) / (2N) apart. 
This inequality follows from the choice of N. Applying (13. 3p and the 
scaling law for capacity yields 

cap (A V:Q ) = a 2 " d cap (aA ViQ ) > a 2 ~ d cX (a(p ) d ~ 2 = cX^ 2 , 

where X = m is the number of sub-cubes Q" of Q in our construction 
that contain at least one point of V . Hence, 

E[cap (A v n Q)\ > c(j) d - 2 E[X}. (3.4) 

The probability that a particular sub-cube Q" contains a point of V is 

l-P(?nQ" = 0) =l-e-^ Q " } , 

by the Poisson nature of the random point process. For any sub-cube 
Q" with centre x, say, 

fi(Q") ~ u(x) ( ~ z/(x)^V choice of N) 

\ N J a 2 

< v{x) £(Q) 2 (f)^ 2 (since a > c/l{Q)) 

= 0(1) (bydEZD). 
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It then follows that 

E[X) = 1 - e " MQ,,) > E * C MQ). I 3 - 5 ) 

Q"CQ Q"CQ 

the last inequality being a consequence of the assumption (11. 5p and 
the fact that the volume of the union of the cubes Q" is some fixed 
fraction of the volume of Q. When combined with ( I3.4p . the estimate 
( 13. 5 p yields the lower bound for E[cap (Ap fl Q)]. □ 



4. Proof of Theorem [Q 
To begin with, we need the following result from Lundh's paper |10j . 

Lemma 4. Let r G S. Then K[W(Ap, r)] is finite if and only if the 
series W(Ap,r) is convergent for almost all random configurations V '. 

Proof. It is clear that E[W(Ap, r)] being finite implies that W(Ap, r) 
is almost surely convergent. The reverse direction is proved by Lundh 
[TO] p. 241] using Kolmogorov's three series theorem. Indeed, it is a 
consequence of this result [SI p. 118] that, in the case of a uniformly 
bounded sequence of non-negative independent random variables, the 
series Ylk-^-k converges almost surely if and only if X]fc^[^fc] is finite. 
As in the proof of Lemma [21 the series W(Ap,r) is split into n series 
W t (Ap,r) = ^fcLi^fc, each of which is almost surely convergent by 
assumption. The random variables Xk in (12. 3p are uniformly bounded. 
It then follows that ^ fc E[X£] = E[^ fc X£] is convergent, that is 
E[W l (Ap, t)] is finite. Summing over i, we find that E[W(Ap, r)] is 
finite as claimed. □ 

Proof of Theorem Q Let us first assume that the finite Poisson bal- 
ayage condition (11. 9p holds for all r in a set T, say, of positive measure 
o~(T) on the boundary of the unit ball and deduce from this that perco- 
lation diffusion occurs. In fact, we will show more - we will show that 
the random archipelago is avoidable with probability one. By Propo- 
sition [TJ we see that E[W / (Ap, r)] is finite for r e T, hence the series 
W(A<p, t) is convergent a.s. for each r G T. For r G T, set 

F T = {V: W(Ap,r) < oo}. 

so that F T has probability 1. We have 



I— [ E[1 F ] dr = E — ^— / 1 F dr 
T) J T [ FtJ [a(T) J T Ft 



o(T) 

from which it follows that j s lF T dr = c r (T) with probability one. 
Equivalently, it is almost surely true that V G F T for a.e. r G T. 
In other words, it is almost surely true that A-p is minimally thin at a 
set of r of positive measure on the sphere S, hence A-p is almost surely 
avoidable by Lemma [TJ 
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Next we prove the reverse implication. For a random configuration 
V, set 

M-p = {r G § : Ap is minimally thin at r}, 

similar to (12. 2p . Suppose that percolation diffusion occurs. Then, with 
positive probability, Ap is minimally thin at each point of some set of 
positive surface measure on the sphere, so that E [J § l Mp {r) dr\ > 0. 
Interchanging the order of integration and expectation, we conclude 
that there is a set T of positive measure on the sphere § such that 
P(r G M-p) > for r G T. By Lemma El W(A v ,t) < oo a.s. for 
r G T. By Lemma H E[W (Ap , r)} is finite for r G T. Finally, it follows 
from Proposition [1] that, for r in the set T of positive measure on the 
sphere §, the Poisson balayage (11.91) is finite. □ 

Proof of CorollaryUl In the case that both <f> and v are radial, the value 
of the Poisson balayage in (11.91) is independent of r 6 S and equals 

[\l - t 2 )0(t) d - 2 ( [ J— ^ da(x)) v(t) dt, 
Jo \Jts \ T ~ x \ J 

where da is surface measure on the sphere £§. Hence (11.91) is equivalent 
to (11 . 10p in the radial setting. □ 



5. Percolation diffusion in space 

The Wiener criterion for minimal thinness of a set A at oo in R d , 
d > 3, is 



, - v -v k 



(see [8], for example) where the cubes {Qk} are obtained by partitioning 
the cube of sidelength 3 J (centre and sides parallel to the coordinate 
axes) into 3 jd cubes of sidelength 3 J_1 and then deleting the central 
cube. Assuming that the radius function <p and the intensity of the 
Poisson process v are roughly constant on each cube Qk (+ version of 
(II. 7p ). the relevant version of Lemma [3] is that, for a cube Qk and any 
point Xk G Q, 

E[cap (A-p n Qk)] ~ (f)(xk) d - 2 fi(Q k ), (5.2) 
and the relevant version of Proposition [1] is 

E[W(A v ,t)]~ [ (^y) u(x)dx. (5.3) 

JR d \B \ \ x \ J 



Since the random archipelago Ap is avoidable in this setting precisely 
when it is minimally thin at the point at infinity, the criterion for 
percolation diffusion is that the integral on the right hand side of (15. 3 p 
be finite. Again this agrees in principle with a criterion for avoidability 
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in the deterministic, regularly located setting Theorem 2] (see also 
[5J Theorem 6]) which corresponds to constant v and 4> radial, namely 

/oo 
r<p(r) d ~ 2 dr < oo. 
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